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Quantification of quantum entanglement in a multiparticle system of two-level atoms 
interacting with a squeezed vacuum state of the radiation field 

Ram Narayan Deb 

Department of Physics, Chandemagore College, Chandemagore, Hooghly, Pin-712136, West Bengal, India 

We quantify multiparticle quantum entanglement in a system of N two-level atoms interacting 
with a squeezed vacuum state of the electromagnetic field. We calculate the amount of quantum 
entanglement present among one hundred such two-level atoms and also show the variation of that 
entanglement with the radiation field parameter. We show the continuous variation of the amount 
of quantum entanglement as we continuously increase the number of atoms from N = 2 to N = 100. 

We also discuss that the multiparticle correlations among the N two-level atoms are made up of all 
possible bipartite correlations among the N atoms. 

PACS numbers: 03.67.Mn, 42.50.Dv, 03.65.Ud 


I. INTRODUCTION 

The quest for the quantification of quantum entangle¬ 
ment in multiparticle systems is an important new arena 
of present day research. Over the past few years there 
has been a growing interest in studyin g quantum en¬ 
tanglement in multiparticle systems [l|-[7|. Peres and 
Horodecki i, 0 made an important step toward the 
understanding of quantum entanglement in the context 
of bipartite states. In Ref. the sufficient and nec¬ 
essary conditions to induce entanglement on two remote 
qubits, by means of their respective linear interaction, 
with a two-mode driving field have been studied. A lot 
of work has also been done to understand the relation¬ 
ship between quantum entanglement and spin squeezing 
[ll[-[2^. In Ref. [2l|, a complete set of generalised spin 
squeezing inequalities for detecting entanglement in an 
ensemble of qudits have been presented. It has also been 
shown, how to detect fc-particle entanglement and bound 
entanglement. In this paper, we quantify quantum en¬ 
tanglement present in a system of N two-level atoms in 
interaction with the squeezed vacuum state of the radi¬ 
ation field. We also discuss that the multiparticle corre¬ 
lations among the N two-level atoms are made up of all 
possible bipartite correlations among the N atoms. We 
also notice that the system shows high value of quantum 
entanglement, when there is no spin squeezing, as defined 
in Ref. [1^, in the system. 

In Ref. [^, we proposed the necessary and sufficient 
condition for the presence of quantum entanglement in 
arbitrary symmetric pure states of N two-level atoms. 
We introduced a parameter, called quantum entangle¬ 
ment parameter, to detect and quantify quantum entan¬ 
glement in such multiatomic systems. In this paper, we 
use that parameter to quantify quantum entanglement 
in a real physical system. The physical system, that we 
take, is a system of N two-level atoms interacting with 
the squeezed vacuum state of the electromagnetic field. 
This system is of considerable interest in Quantum Op¬ 
tics. We quantify the amount of quantum entanglement 


among the atoms of this system for the number of atoms, 
varying from N = 2 to N = 100. We also show that the 
system of atoms shows high quantum entanglement when 
there is no spin squeezing, as defined in Ref. [j^ . in the 
system. We emphasize the fact that, the quantum entan¬ 
glement in such multiatomic systems can be quantified 
experimentally. 

In section II, we present the construction of the quan¬ 
tum entanglement parameter. In section III, we present 
our study on quantum entanglement of N two-level atoms 
interacting with the squeezed vacuum state of the radi¬ 
ation field. In section IV, we present the summary and 
conclusion. 


II. CONSTRUCTION OF QUANTUM 
ENTANGLEMENT PARAMETER 

An atom has many electronic energy levels, but when 
it interacts with an external monochromatic electromag¬ 
netic field, the atom makes a transition from one of its 
energy level to the other. In this case, we mainly con¬ 
centrate on those two energy levels and hence the atom 
is called as a two-level atom. 

We consider a system of N such two-level atoms. Now, 
if among the assembly of N such two-level atoms, the n- 
th atom has the upper and lower energy levels, denoted 
as \un) and \ln), respectively, then, we can construct the 
pseudo-spin operators (with h=l), 



= (1/2)(|m„)(/„| -1- ), 

( 1 ) 

Juy 

= (-z/2)(|m„)(;„| - |;„)(u„|), 

( 2 ) 

Jriz 

= {l/2){\Un){Un\-\ln){ln\), 

( 3 ) 

such that, 


Jny\ — 

( 4 ) 


and two more relations with cyclic changes in x, y and 
z. For the entire system of N two-level atoms, we have 
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collective pseudo-spin operators, 

N N N 

Jx = ^ ^ Jia ;) Jy = ^ ^ Jiy ) >^ 2 : = ^ ^ Jiz ■ ( 5 ) 

i=l i—1 i—1 

The individual atomic operators satisfy 

[ii„, ^2,,] = 0, [jl^,Jly] = [J2^, J2y] = iJ2^, ■■■{&) 


As a consequence of these commutation relations, the 
collective pseudo-spin operators Jx-, Jy and Jz satisfy. 


\Jx, Jy\ — iJz 

(7) 

and two more relations with cyclic changes 

in X, y and z. 

The simultaneous eigenvectors of J^ = 

jj + jy + Jz 

and Jz are denoted as \j, m) where 


=j{j + l)|j.w) 

( 8 ) 

and 


II 

(9) 


and for an atomic squeezed state, also called squeezed 
spin state (SSS) (^ . 

i N 

AJx>^ or, AJy,^<l = -. (14) 

It is to be mentioned here that 

N 

Ajy,,AJy, > — , (15) 

which is the Heisenberg’s uncertainty principle. 

A collective state vector |a) for a system of two atoms 
is said to be quantum mechanically entangled if |a) can¬ 
not be expressed as a direct product of the two individual 
atomic state vectors i.e. 

\a) ^ \ai) ® \a2): (16) 

where |ai) and J^) are the state vectors of the two in¬ 
dividual atoms [23|. 

We now proceed to construct the quantum entangle¬ 
ment parameter for a system of N two-level atoms. 

An arbitrary symmetric pure state for a system of N 
two-level atoms in the {toi, m 2 , TO 3 , ....WAr} representa¬ 
tion is given as 


The quantum number j is related to the number of atoms 
N as j = N/2 and m = —j,—j + l,----(j — l),j- The 
collective quantum state vector for a system of N two- 
level atoms can be expressed as a linear superposition of 
\j,m) as 
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IV'i) = 51 Cm|j,m). (10) 

m=-j 

Now, to know whether an arbitrary quantum state vector 
lipj) of the system is an atomic coherent state or atomic 
squeezed state, we calculate the mean pseudo-spin vector 

{3) = {Jx)i+{Jy)j + {J.)k, (11) 

where i, j and k are the unit vectors along positive x, 
y and z axes respectively, and the average values in the 
above expression are to be calculated over the state {ipj). 
The mean pseudo-spin vector (J) points in an arbitrary 
direction in space. Therefore, we conventionally rotate 
the coordinate system {x,y,z} to {x',y',z'}, such that 
(J) points along the z' axis and calculate the quantum 
fluctuations in J^' and Jy' for the state I'i/’j)- These quan¬ 
tum fluctuations are 


AJx',y' — ii^jlJx'{'^j\Jx',y'\'4’j) ■ (12) 


Now for an atomic coherent state, also called coherent 
spin state (CSS) in literature, 22|, the above quantum 
fluctuations are equal and they are, 
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J_ 

2 


N 

T’ 


Ivk) 


111 i\ 

rrr"" 2 / 


1 

1 1 

1\ , 

1 1 1 

iV 

2’ 

~2’2’' 

- 2 / + 

2’2’2’'' 

2 / 


= Gi 


1 1 
2 ’ 2 ’ 


Go 




+ 


+ 


G 3 

V ^2 


111 1 \ 

r~rr"" 2 / 


1 

1 

1 

1\ 

\ J_ 

1 

1 

1 

iV 

~ 2’ 

2’ 

~2’' 

■'2/ 


2 

2 ’" 

■■ ~ 2’ 

2 /. 


■Gn+1 


1 1 1 
2 ’~ 2 ’~ 2 ’ 



(17) 


where Gi, G 2 , 
as 


..., Gat+i are constants and ^Cr is given 


N 


Cr 


Nl 

r\{N — r)! 


(18) 


Now, the mean pseudo-spin vector (J) for the above 
state points in an arbitrary direction in space. Assuming 
that it lies in the first octant of the coordinate system, we 
perform a rotation of the coordinate system from {x, y, z} 
to {x',y',z'}, such that (J) points along the z' axis. The 
operators {Jx', Jy', Jz'} in the rotated frame {x',y',z'} 
are related to {Jx, Jy, Jz} in the unrotated frame {x, y, z} 
as 


Jx' = Jx cos 0 cos 0-I-jy cos 0 sin (/) — jj sin 0 (19) 
Jy! = —JxS\n(j) + JyCOSCj) (20) 

Jz' = JxSi'!i9cos(j) + JySinOsia(j) + JzCOs9, (21) 


(13) 
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where, 


Using Eqs. m, we also have 


cos 9 


COS(j) 



( 22 ) 

(23) 


We can observe that for the above choice of cos 9 and 
coscj), we have {Jx') = 0, {Jy>) = 0, and the mean spin 
vector points along the z' axis. 

Now, using Eqs. (fTO|) and (1^ the quantum fluctua¬ 
tions A Jj, and A J^, for this system are, 


-f 

X 


AJJ, = (J2,) - (J,,)2 

AJx^ cos^ 9 cos^ (f) + A Jy"^ cos^ 9 sin^ ( 

AJz^ sin^ 9 + ^{JxJy + JyJx) ~ 2( Jx 

cos^ 9 sin (j) cos </> — ({JxJz + JzJx) — 


X sin 9 cos 9 cos </) — [^{JyJz + JzJy) 
X sin 9 cos 9 sin (j) 



{JxJy + JyJx) 2(Jx)(Jy) — 'y ) “f 

i=l 

N N 

- 2(4)(4)) +2^^((4j,J-(4)(4)) (32) 




and 


and 


A4 = (4) - {Jy4 

= AJx^ sin^ (j) + AJy^ COS^ <() — ( {JxJy + JyJx) 


lx)\Jy) I Sm(pCOS(p. 

Now using Eqs. ®, we have 


and 


N 


N N 


4 - X! 4» + X! X! ’ 
2=1 2=1 1 = 1 

N N N 

4 " = T.Jl + Y.T.J-A, 

2=1 2 = 1 1=1 

N N N 

4 = E 4 + EE 44 


2=1 2=1 l = j 


N 


N N 


2 = 1 




(25) 

(26) 

(27) 

(28) 


(Jx)^ = E(4)^ + EE(4)(4) (29) 

2 = 1 2=1 1 = 1 

N N N 

(Jy)" = E(4)'+EE(4)(4) (30) 

2=1 2=1 i = l 

N N N 

(j4 = E(4)^ + EE(4)(4)- (31) 


{JxJz + JzJx) — ^{Jx){Jz) = + JizJix) 

2=1 

N N 

- 2(J2j(4))+2^^((44)-(4)aj) (33) 


2=1 i = l 


and 


{JyJz + JzJy) ‘^{Jy){Jz) — ^('^4'^4 + ^iz^iy) 

2=1 

N N 

- 2(4)(4)) +2EE((44)-(4)(4))-(34) 
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Now, using Eqs. ([26|) - (|34|) in Eq. (|24|) . we obtain 


AJi = 


N N 


N N N N 


SpJ 


+ 


+ 


+ 


+ 


El: (4) (4)} 


2=1 1 = 1 
l^i 


COS^ 0 COS^ (j) + 
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E(4> 


N N N 

E E ) “ {E ^'^4 ^ 
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■ ^ ' 2^1 2=1 1 = 1 

N 
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X sin 9 + 


N N 
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2EE{(EE)-(4)(E)} 




AT 


cos^ 9 sin (j) cos (f) 


E{(E4 + 44)-2(4)(4)} 


2=1 
N N 
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+ 


2EE{(E4)-(4)(4)} 

1 ;' 

N 

E! I ("^4 "^2 + "^2 Jiy ) — ‘^{Jiy ) {Jit ) I 

i=l 
N N 

2EE{(44)-(4)(4)} 


sin 9 cos 9 cos (j) 


2=1 1 = 1 
l^i 


sin^cos^sint/). 

(35) 


The above expression can be rearranged and written as 


N 




AJ~ cos 9 cos (j) + AJ~ cos 9 sin (j) 


X sin 9 cos 9 sin (f) 


N N 

EE 




X cos^ 6 »cos^^+ - {Ji^){Ji^)^ cos^ 9s\T? 4> 

+ ((44) - (4)(4)) sin2 0 + 2((44) 

- {Jix){Jiy)) COS^ 0 sin COS - 2((4- 2 / 2 ) 

- (4 ) {Jix )) sin 9 COS 6* COS (/> - 2 ((4 ) - (4) (>^4)) 


X sin 9 cos 9 sin (j) 


(36) 


If, we now compare the term under the single summation 
symbol with Eq. ((Ml) , we find that it is A J?,. Therefore, 
Eq. (1551) can be expressed as. 


N 


AJl, = E A 4 , + CORRX, 


(37) 


i=l 


where CORRX is the term with the double summation 
symbol in Eq. (l36l) . We notice that CORRX is solely 
made up of the quantum correlations among the N two 
level atoms. Thus, the quantum fiuctuation A Jj, for the 
composite state of N two-level atoms has been expressed 
as the algebraic sum of the quantum fluctuations AJf ^ 
of the N individual atoms and the quantum correlation 
term CORRX among the N atoms. We now present a 
similar expression for AJ^,. Using Eqs. (1551) . (1571) . (1551) . 
(1551) . and (1551) in Eq. (1551) . we obtain 


AJ^, = 

y 


N N 


N N N N 

E(4)+EE<44)-{E<4)" 


2=1 


EE<4)(4)} 


2=1 i = l 


N N 


2=1 1 = 1 


Sin 
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N N 
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N N 


-f 


2EE{(E4)-(4)(4)} 




sin (/) COS (/>. (38) 


+ ^4 ^ + ((4 

Jiy JiyJix) J'{Jix){Jiy)'^ 

> 

II 

M* 

X cos^ 9 sin (j) cos (j) — 

f(44 + 44)-2(4)(4)) 

2=1 _ 


The above expression can be rearranged as 
Ajf sin^ 6 + A jf cos^ i 


X sin0cos0cos0 - ((44 +“ 2(4)(4)) 


{^{^ix^iy Jiy^ix) ‘^0ix){Jiy)^ sini^COSl^ 
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N N 

EE 


JiJiJ - (E)(-^Zx)) sin 2 </> 


+ ((44) - (4) (4)) cos4 + 2((44) 

- ( 4 )( 4 )) sin</'COS 0 . 


(39) 


Comparing the term under the single summation symbol 
in Eq. (l3^ with the right hand side of Eq. (|2^ . we 
notice that it is AJf ,. Thus, Eq. (l3^ can be written as 

N 

A A J2 ^ + CORKY, (40) 

i=l 

where, CORKY represents the term with the double 
summation symbol in Eq. (15^ . We notice that CORKY 
is solely made up of the quantum correlations among the 
N two-level atoms. Thus, the quantum fluctuation A J^, 
for the composite state of N two-level atoms has been 
expressed as the algebraic sum of the quantum fluctua¬ 
tions A 4 of the N individual constituent atoms and 

^y' 

the quantum correlation term among them. 

Now if the quantum state is an unentangled state, 
that is 


l^*) = l^'i) ® |4'2) ® l^-a).Ol^-jv), (41) 

where I'I'i), |'p 2 ), l^'iv) are the individual atomic state 
vectors, then, we have for i ^ I, 

(44) = (4)(4), (44) = (4)(4)> 

{JizJlz) — 0iz)0iz)'> 0ix'^ly) — {•Jix)0iy)'> 

{^ixJlz) — {Jix){Jlz)i (JiyJlz) — 0iy){Jlz)- (^ 2 ) 

In this case the terms under the double summation sym¬ 
bols in Eqs. (|36l) and (|39|) are zero, implying that 

CORRX = CORKY = 0. (43) 

Thus, for an unentangled state we have CORRX = 
CORKY = 0. 

Now, we can calculate and find that the quantum fluc¬ 
tuations and AJ^^, ^ of the individual constituent 

atoms have value j, that is 



II 

(M 

<1 

II 

i> 

sy 

to 

1 

“ 4’ 

(44) 

AJi 

Xy ! 

II 

> 

to 

II 

i> 

to 

1 

“ 4‘ 

(45) 


Now, it has been discussed that [l^, N/A is the quantum 
fluctuation in J^' and Jy/ for a system of N two-level 
atoms in a coherent spin state (CSS), which is a com¬ 
pletely separable state, that is, an unentangled state. So, 
we can write Eqs. (H51) and (H71) as 

AJy,,^ - AJy,.Xn-ent. = CORRX, (48) 

AJy,^ - AJy,Xn-ent. = CORKY. (49) 


Therefore, CORRX and CORKY are the measures of 
the deviations of the quantum fluctuations in J^i and Jyi 
respectively, of the state |4>) from those of an unentangled 
state. 

If for some quantum state of the N two-level atoms 
we have CORRX < 0, then AJX < A/4, and we say 
that the corresponding quantum state has spin squeez¬ 
ing in the x' quadrature. In that case we must have 
CORKY > 0, as to restore the Heisenberg’s uncertainty 
principle (Eq. (I15|) b Similarly, when CORKY < 0, mak¬ 
ing CORRX > 0, we say that there is spin squeezing in 
the y' quadrature. 

The non-zero value of either CORRX or CORKY or 
both implies the presence of quantum correlations among 
the atoms, because in that case the conditions in Eq. 
(l42l) are not all satisfied and the terms under the double 
summation symbols in Eqs. (IMl) and (123) are non-zero. 
Now, it may happen that for some quantum state of N 
two-level atoms, both CORRX and CORKY are greater 
than zero, implying no spin squeezing at all, but since 
CORRX and CORKY are non-zero, there is quantum 
correlations among the atoms and we have entanglement. 
Thus, the presence of entanglement cannot be detected 
always by spin squeezing and we show this for a real 
physical system in the next section. 

We notice from Eqs. (I5S|) - (1271) . (1221) . and (03, that 
CORRX and CORKY are the extracts from the com¬ 
posite quantum fluctuations AJ^, and AJy, of the N 
two-level atoms, representing only the quantum correla¬ 
tions among the N atoms. We also observe from the 
above mentioned equations that CORRX and CORKY 
are made up of all possible bipartite quantum correlations 
among the N atoms. Therefore, we use these correlation 
terms, CORRX and CORKY, to construct the multi¬ 
particle quantum entanglement parameter. Now, since, 
CORRX and CORKY may be positive or negative, we 
construct the multiparticle quantum entanglement pa¬ 
rameter E as 


Therefore, using Eqs. (1221) and (1251) in Eqs. (1271) and 
(03, we obtain 


AJX = —+CORRX, (46) 

4 

AJy,'^ = CORKY. (47) 


E = - 
2 1 


{CORRXf Y {CORKYf 


(50) 


Thus, E is the mean squared deviation of the quantum 
fluctuations in J^' and Jy' of an arbitrary state from 
those of an unentangled state. Now, using Eqs. (1251) and 
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621 ), we can express E as, 

1 


E = 


2 L 




(51) 


Now, if and are the spectrosco pic squeezing pa¬ 
rameters used in Ramsey spectroscopy [2^, given as 


= ^AJ-' . A, = ^AJ„, 


then. 


E = 


|(J)| |(J)| 

C|J(J)IV?LKJ)I" N 


y' ■> 


(52) 


2j V 2j 2 
enji^A/aji^A N 


2 j 


2 j 


Ar 2 


(53) 


after the operator in Eq. ((Ml) , we get 

l^'m) = \j,m'){j,m'\e~"^-^«\j,m) 

m'^-j 

= Y ^m'm(V 2 )|j,w') 

m'^-j 

= Ai Y e^'"'^m'm(V 2 )|j,TO'), ( 56 ) 


m'^-j 


where d^/m('^/2) is the reduced Wigner d-matrix element 
[i^ given as, 

dL'AP) = {E'm'\Ry{P)\j,m) = (j, |j, m) 


= (- 1 )™' + w)!(j - ”^)! 0 ' + ™')Kj - "i')! 


Thus, the quantum entanglement parameter E is con¬ 
nected to the experimentally measurable quantities. 

Now, if for a quantum state of N two-level atoms, 
E = 0, then, both CORRX and CORRY are zero, that 
is all the conditions in Eqs. (|42l) are satisfied and the cor¬ 
responding quantum state is completely separable, and 
we have no entanglement. If E > 0, then either CORRX 
or CORRY or both are non zero, and the correspond¬ 
ing quantum state is an entangled state. Now, the nu¬ 
merical value of E can be a measure of the amount of 
quantum entanglement present among the N two-level 
atoms. Thus, by measuring E with the help of Eqs. m 
or (IM|) for a quantum state of N two-level atoms we can 
measure the amount of quantum entanglement present in 
that system. In the following section we apply this idea 
to measure the amount of quantum entanglement present 
in a real physical system, which is of considerable interest 
in Quantum Optics. 


III. QUANTUM ENTANGLEMENT OF A 
SYSTEM OF N TWO-LEVEL ATOMS 
INTERACTING WITH THE SQUEEZED 
VACUUM STATE OF THE RADIATION FIELD 

We consider a system of N two-level atoms in interac¬ 
tion with the squeezed vacuum state of the electromag¬ 
netic field. The composite quantum state vector for such 
system is given as [i^, [ 2 ^ . 

I^*™) = (54) 

where is the normalization constant, ^ is the radia¬ 
tion held parameter and |j, m) is the Wigner state. In¬ 
serting the identity operator 

Y = / (55) 

m'^-j 


^ (-I)'=(C0S f )2t-2fc-™'+"*(sin Ifk+m'-TU 

k\{j — m' — k)\{j + m — k)l(m' — m + k)l 


■ (57) 


The normalization constant Am is given by 




(58) 


Using the addition theorem and symmetry properties of 
Wigner rotation matrices [i^ , we can obtain 


A-2 = di^mirn = A, 


(59) 


where = A, is the reduced Wigner d-matrix 

element for imaginary angle 2i^, and is given as 


A = (j +m)l{j - m)!' 


(cosh^) 2 J(tanh.f) 


2k 


^ (fc!)2(j — TO — k)\{j Ym — k)\' 

(60) 

In Refs. [ 2 ^ and (2^, the spin squeezing properties of 
the state |4>m) have been studied. But, quantum entan¬ 
glement does not have any relationship with spin squeez¬ 
ing, which we show below. In this paper, we study the 
quantum entanglement property of this state. 

Now, to calculate and quantify the amount of quan¬ 
tum entanglement among the N two-level atoms in this 
state, we need to calculate the quantum entanglement 
parameter E for this state. This in turn, requires the 
calculation of the quantum fluctuations AJ^, and AJ^, 
for this state in a coordinate frame in which the mean 
pseudo-spin vector (J) points along the z'-axis. There¬ 
fore, we now calculate AJj, and AJ^, for the state |4>m). 
Now, |4'm) is an eigenvector of a non-Hermitian operator 
A with eigenvalue to, 25 |, [ 2 ^ as 


A|4'm) = (Jxcosh^ -I- iJy sinh^)|^'m) = m\^m)- (61) 
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The dual of the above is written as 


or, 


(vE'mIA'I' = (TmKJx cosh^ - zJj^sinh^) = {^m\m. (62) 

Taking the scalar product of both sides of Eq. (ED) by 
('kmI and equating the real and imaginary parts we get, 


^ ^ A ‘2iTn '' 

m\JxJz JzJx\^rn) — rn\Jz\^rn) 

cosh 4 


sinh^ 

cosh^ 


(^m \Jx l^m) • 


(75) 


^ jji 

^ (63) 

cosh^ 

and 

=0. (64) 

Taking the scalar product of Eq. (EB) with it’s dual i.e. 
Eq. (1^ viz. 

(4'm|A'l'A|4'm) = TO^, (65) 


we obtain 


cosh^ + sinh^ |4'm) 

+ *sinh^cosh^('I'm|[Tx, JyJl'I'm) = (66) 


or, 

cosh^ f( 4 'm| J^| 4 'm) + sinh^ ^{^rn\Jy\'^m) 

- sinh^coshf(4'm|i2|4'm) = (67) 

Operating with A on both sides of Eq. EB and then 
taking the scalar product with ('km| we get 

(^m|A2|4-m) =m2. (68) 

More explicitly it is 


Using Eq. EB we have 


{'^Tn\JxJz~^ JzJx\'^ni) — 


2m 


— m- 


cosh^ 
sinh^ 


(4'm|Jz|4'm) 


(76) 


cosh ^ 

Subtracting Eq. (1751) from Eq. (17^ we have 

(^m|['72;, Tcjl^km) COSh^ 

+ isinh^(4'm| iyiz + izJy|4'm) = 0 (77) 


or, 

(ikm I <72/1 ikm) COSh^ T sinll ^(\km | *7y T JzJyY^ rn) — 0- 

Using Eq. (IMl) we obtain 

('I'm|"7y>72 + Jy I'f'm) = 0. (78) 

To calculate ('kml^zl'I'm) we proceed as below. 

= JzAm E 

m'^-j 

= Am (79) 


cosh^^(4'm|ixl^m) - sinh^^(Tm|iy |4'm) 

+ *sinh,fcosh^(4'm|ixiy + ^yTa,|4'm) = (69) 

Equating the real and imaginary parts we obtain 

cosh^^(4'm|iJ|4'm) - smh^5(4'm|iy|4'm) =TO^ (70) 
and 

{'^m\JxJy + <7y'7a;|'I'm) = 0. (71) 

We now operate both sides of Eq. (ICT) by Jz from left 
and then take the scalar product with ('I'm|, and obtain 

m\Jzi Jx cosh ^ iJy sinh^)|'km) — rn('km|'7z|'km)- 

(72) 

Taking complex conjugate of both sides, we obtain 

m\{Jx cosh ^ ijy sinh ^),75; |lkm) — rn('I/m|'7z|'km)- 

(73) 

Adding Eqs. (17^ and (1751) . we get 

COSh^(\km|'7zTx T JxJzl'^m} 

+ isinh^(4'm|[iz, Jy]\'^m) = 2m(4'm| |4'm) (74) 


Taking scalar product by we obtain 


m' = —j m" — —j 

m' = -j 

^ m' = -j 

= \Am\^~di,mm)- (80) 

Using Eq. EB we obtain 


(4'm|iz|4'm) = 


1 dA 


2A 

With this expression of {Jz), Eq. (176)) becomes 


(lkm|'7a:'7z T JzJx\^m) — 


m 1 dA 


sinh,^ 


cosh 5 A cosh^^ 


(81) 


• (82) 









In the same manner it can be shown that 


Putting /3 = 2i^, we obtain 


■ (83) 

Adding Eqs. (p71) and (1701) and using Eq. (ISTl) we get 
^2 1 , dA 


= 4j-(j + 1)A-4—^-2coth2e—. (92) 

cosh ^ d? 




m 




(84) 


Using Eq. (IM|) this can be expressed as 


de 


= 4j(j + 1)A - 4 


m?A „ cosh 2^ 


-r- 


cosh^^ cosh^^ 


(93) 


Subtracting Eq. (1701) from Eq. (1071) and using Eq. (10T1) 
we obtain 


( 4 'm|iy| 4 'm) = ^coth^^. 


The quantity ^ is expressed in suitable form. Differ¬ 
entiating once the expression of A as given in Eq. (1001) . 
we obtain 


.dA 


(85) 


Thus, we have obtained all the necessary averages of the 
pseudo-spin operators over I'l'^) to calculate AJ^A and 

Ajy! . 

Now, from Eqs. ( 1001 ) . (| 00 |) and it is evident that 
the mean pseudo-spin vector (J) is not along the z-axis 
and lies in the z — x plane making an angle say 0 i with 
the z-axis. Therefore, we align the vector (J) along the 
z'-axis by performing a rotation as below. 


rfA 

— = tanh^r 

(86) 

Jx' 

= Jx cos 01 — Jz sin 01 

(94) 

di 

Jy' 

= Jy 

(95) 

with 


Jz' 

= Jx sin 01 -k Jz cos 01 

(96) 

r = 2j A -k 2 ^ „ 

cosh^e 

(87) 

with 

tan 01 = IE, 

0^) 


where ry is given as 



(97) 


77 = (cosh^)^'^ (j-I- 77 i)!(j — m)! 

(tanh^)^'^ 


E 


k\{k -I- l)!(j -I- m — 1 — k)\{j — m — 1 — k)\' 


( 88 ) 


This rotation makes {Jx') = 0 and since {Jy) is already 
zero as given in Eq. (1001) , the vector (J) is now along the 
z'-axis. We now observe the variances in Jx' and Jy'. As 
Jy! = Jy, hence the variance in Jy' is the same as Jy i.e. 


l2 A 

To express in suitable form we use the differen¬ 
tial equation satisfied by the rotation matrix element 
Z3^,j^(a,/3, 7 ) familiar from the quantum mechanics of 
a symmetric top [ 2 ^ . The rotation matrix element 
D^/^(a,/ 3 , 7 ) is defined as 


= \J{Jp) - {Jy-) 


— Ajy. 


(98) 




(89) 


As (Jy) = 0 by Eq. (1001) . therefore. 


where R{a, /3, 7 ) is the rotation matrix and a, /3 ,7 are 
the Euler angles. Now, the differential equation satisfied 

by is 


19 — 2 mm'cos ^ 

;Sinp- 


AJy" = 01) = ^ cothe^ = AJy,^ (99) 


where we have used Eq. (1501) . Using Eq. (1501) . we get 

r 


sin^ /3 


A T^, = 

'^y' 4A- 


( 100 ) 


sin df3 djd 
X =j(j +l)^Lm'(«>/?> 7)- (90) 

Now, by putting m' = m and a = 7 = 0 we get 
d2 


On the other hand since {Jx') = 0, the square of the 
variance in Jx' is given as 


dp- 


-dlnmiP) = -j'(j + l)c?mm(^) + W SeC (/3/2)d^„ (/3) 


- cot/3—d^„„(/3). 


(91) 


AJj, = 01 '). 

Therefore, using Eq. ([M]). we get 

Aj'^i = (Jj) cos^ 9 i + sin^ 61 
- {JxJz + JzJx) sin 0i cos 01 . 


( 101 ) 


( 102 ) 
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Using Eqs. (|82]), ([83]), (|84|), ([ 86 ]), ([93]) and ([97]) in Eq. 
(I102|) . we obtain 


AJ2 = 


m 


cosh ^ 


tanh^ £ / r \" 

+ ^(a) 


^tanh^ ^ 
V 4 


(S' 


j(j + l)m^ r 

cosh^ £ 4 cosh^ £ A 


cosh^ £ 


(103) 


It is to be noted that since ■£ is symmetric under m —>■ 
—TO, which can be verified from the Eqs. ([60]), ([87]) and 
([ 88 ]), the quantity AJ^, and AJ^, are also symmetric 
under to —> —to. For to = ±j we have 


AJS = A4 = 


(104) 


The value of |(J)| for the state I'h^) is obtained using 
Eqs. (1551) . (IMl) . (1^ and (1551) as 


KJ)I = 


Using Eqs. (llOOl) and (|103l) in Eq. (1511) . we can calcu¬ 
late the value of the quantum entanglement parameter E, 
numerically, for various values of the radiation field pa¬ 
rameter £, and for various numbers N of two-level atoms 
in the system. The numerical value of E can be taken 
as the amount of quantum entanglement present among 
the atoms. If E is zero, then there is no quantum entan¬ 
glement among the atoms and the atoms are completely 
uncorrelated. 

Using Eq. (11051) in Eq. (l53l) . and experimentally mea¬ 
suring the spectroscopic squeezing parameters £/{^ and 
^Ry j can calculate the quantum entanglement param¬ 
eter E from Eq. (ES]). 

We, now, show in Fig. 1, the variation of the quan¬ 
tum entanglement parameter E for one hundred two-level 
atoms in interaction with the squeezed vacuum state of 
the radiation field. We plot graphs for various values of 
TO. The non-zero value of E indicates the presence of 
quantum entanglement and the numerical value of E can 
be taken as the amount of quantum entanglement present 
among the atoms of the system. 

We observe from Fig. 1 that the quantum entangle¬ 
ment parameter E starts from a very high value when the 
radiation field parameter £ = 0 , and sharply decreases, as 
£ increases. We can find that E is very close to 0, when 
£ is close to 3. We also observe that, for a fixed value of 
£, the amount of quantum entanglement increases with 
decrease in to values. We note that, E is always zero for 
TO = 50, which is the maximum value of to for N = 100. 


^( J ,)2 + ( J ,)2 + ( J ,)2 

■ m? tanh^ ^ ^ 

- cosh^ £ 4 V A / 


2i 


1/2 



FIG. 1: Variation of the quantum entanglement parameter 
U as a function of the radiation field parameter ^ for N = 
100. E and ^ are plotted on the vertical and horizontal axes 
respectively. 


TABLE I: Values of E with ^ for different values of m, with 
N = 100 


m C 

E 

£ 

E 

£ 

E 

£ 

E 

10 0 

1440000 

0.1 

22357.14 

0.2 

4914.34 

3 

0.0142 

20 0 

1102500 

0.1 

19162.98 

0.2 

4189.47 

3 

0.0109 

30 0 

640000 

0.1 

13943.33 

0.2 

3013.98 

3 

0.0063 

40 0 

202500 

0.1 

6965.37 

0.2 

1470.5 

3 

0.002 

50 0 

0 

0.1 

0 

0.2 

0 

3 

0 


In Table I, we present the values of E with £ for different 
values of to, when N = 100. 

The numerical values of E, given in Table I, can be 
taken as the measure of quantum entanglement for the 
corresponding values of £, for hundred atoms. 

In Fig. 2 we graphically show the variation of the 
quantum entanglement parameter E with the number of 
atoms varying from A = 2 to A = 100 for fixed values 
of £. Here, we have taken to = 1 . The two curves cor¬ 
respond to two fixed values of £, which are £ = 0.8 and 
£ = 1. We observe that E increases with A, and also the 
values of E are higher at £ = 0.8 than those at £ = 1. We 
can also see from Eqs. (I100|) . (11031) . and (11051) . and from 
Eqs. (IKTl) to ([53]) . that the amount of quantum entan¬ 
glement present in the system is independent of the sign 
of TO, that is, E is invariant under the transformation 

TO - > —TO. 

In Fig. 3 we show the variation of A.Jx'^ and AJyi'^ 
with the radiation field parameter. In the figure, these 
quantities have been shown as Var — Jx and Var — Jy 
respectively. We observe that both the quantities AJ^,/^ 
and AJyi'^ have values greater than 100 in the range from 
£ = 0 to £ = 0.03. Thus, in this range of £, the quantum 
state shows no spin squeezing at all, whereas, we 

observe from Fig. 1 that in this range of £, the quan- 
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FIG. 2: Variation of the quantum entanglement parameter E 
with the number of atoms N for two fixed values of which 
are ^ = 0.8 and ^ = 1. Here we have taken m — 1. E and N 
are plotted on the vertical and horizontal axes respectively. 



FIG. 3: Variation of and AJyi^ with radiation field 

parameter ^ for m = 40. and AJyi^ have been shown 

as Var — Jx and Var — Jy respectively in the figure. 


turn entanglement parameter E has a very high value, 
implying the presence of quantum entanglement among 
the two level-atoms. Thus, presence of quantum entan¬ 
glement among the atoms does not ensure spin squeezing, 
as defined in Ref. [l^, of the system. 


IV. SUMMARY AND CONCLUSION 

We quantify the amount of quantum entanglement 
present in a system of N two-level atoms in interaction 
with the squeezed vacuum state of the electromagnetic 
field. The method, that we use to quantify the amount 
of quantum entanglement present in such system, is to 
calculate the numerical value of the mean squared devi¬ 


ations of the quantum fluctuations in x' and y' quadra¬ 
tures of the considered quantum state from the corre¬ 
sponding quantum fluctuations of an atomic coherent 
state or coherent spin state (CSS). We use this method, 
because, we have expressed the composite quantum fluc¬ 
tuations in x' and y' quadratures of N two-level atoms as 
an algebraic sum of the corresponding quantum fluctua¬ 
tions of the N individual constituent atoms and the quan¬ 
tum correlation terms among them. We have found that 
these multiparticle quantum correlation terms CORRX 
and CORRY are equal to the deviations of the quantum 
fluctuations in x' and y' quadratures respectively, of the 
corresponding composite quantum state from those of an 
atomic coherent state or coherent ™in state (CSS), which 
is a completely separable state 0, [i^. We have also 
seen that the multiparticle quantum correlation terms 
CORRX and CORRY are made up of all possible bi¬ 
partite (two-atoms) quantum correlations among the N 
atoms. Now, since these multiparticle correlation terms 
may be positive or negative, we take the mean of the 
squares of CORRX and CORRY as a measure of the 
amount of multiparticle quantum entanglement present 
in the system of N two-level atoms. This mean squared 
value of CORRX and CORRY is equal to the mean 
squared value of the deviations of the quantum fluctua¬ 
tions in x' and y’ quadratures of the considered quantum 
state from the corresponding quantum fluctuations of an 
atomic coherent state. We, thus, quantify the amount of 
quantum entanglement present in the system of N two- 
level atoms. 

We graphically show the variation of the multiparticle 
quantum entanglement with the radiation field param¬ 
eter for one hundred such two-level atoms, for different 
values of m. We note that the amount of quantum entan¬ 
glement for a fixed value of the radiation field parameter, 
increases with the decrease in the value of m. We also 
note that, the amount of quantum entanglement present 
in the system remains invariant under the transformation 
m —> —m. We show graphically, the continuous varia¬ 
tion of the amount of quantum entanglement present in 
the system, as we increase continuously the number of 
atoms in the system from N = 2 to N = 100, for two 
fixed values of the radiation field parameter. We note 
that, the amount of quantum entanglement for any num¬ 
ber of atoms N is higher for lesser value of the radiation 
field parameter. We also note that, the amount of quan¬ 
tum entanglement in the system of atoms increases, as 
we increase the number of atoms in the assembly. 

We also show that in certain range of the radiation 
field parameter, the system shows no spin squeezing, as 
defined in Ref. [l^, but shows high quantum entangle¬ 
ment. 

Our method can also be applied to the system stud¬ 
ied by Felicetti et ah, [s^, where they propose a super¬ 
conducting circuit architecture for multipartite entangle¬ 
ment generation. It can also be applied to the system 
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studied in Ref. 31|, in which they propose a scheme 
to controllably entangle the internal states of two atoms 
trapped in a high-finesse optical cavity by employing 
quantum reservoir engineering. Our method can also be 
applied to the system used in Ref. [s^ , where they stud¬ 
ied the dynamics of the entanglement of formation in a 
quantum-measurement model, consisted of a four-level 
atom harmonically bound in a three-dimensional trap. 

We hope that our method of quantification of the 
amount of quantum entanglement in multiparticle sys¬ 
tems may produce some new insight into the subject. 
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